Particle-number scaling of the phase sensitivity in realistic Bayesian twin-mode 

Heisenberg-limited interferometry 
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We investigate the scaling of the phase sensitivity of a nonideal Heisenberg-limited interferometer 
with the particle number A*", in the case of the Bayesian detection procedure proposed by Holland 
and Burnett [Phys. Rev. Lett. 71, 1355 (1993)] for twin boson input modes. Using Monte Carlo 
simulations for up to 10 000 bosons, we show that the phase error of a nonideal interferometer scales 
with the Heisenberg limit if the losses are of the order of A^~^. Greater losses degrade the scaling 
which is then in N~^^'^ , like the shot-noise limit, yet the sensitivity stays sub-shot-noise as long 
as photon correlations are present. These results give the actual limits of Bayesian detection for 
twin-mode interferometry and prove that it is an experimentally feasible scheme, contrary to what 
is implied by the coincidence-detection analysis of Kim et al. [Phys. Rev. A 60, 708 (1999)]. 



I. INTRODUCTION 

Heisenberg-limited quantum interferometry (HLI) is 
of interest for ultra-precise measurement of phase and 
energy differences Possible applications include 

the enhancement of atomic frequency standards and of 
gravitational-wave interferometers. Although optical in- 
terferometry often comes to mind first, interferometry at 
the Heisenberg limit concerns any bosonic quantum wave, 
and even fermionic ones . The absolute minimum error 
on the phase difference (j)- between the two interferome- 
ter arms, for TV detected bosons, is the Heisenberg limit 
(HL), 



A(j)_ - — , 
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to be compared to the shot-noise limit of the beam split- 
ter (SNL), 
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These limits arise from the quantum interference of the 
two input modes, a and 6, on the first beam splitter 
(boson-mode splitter) of the interferometer (Fig.Q). This 




V2 



(a + ib) 



V2 



b) 



interference defines the quantum noise AA^_ = Na — Ni, 
of the photon number difference at the output of the 
beam splitter (i.e. inside the interferometer). The quan- 
tum noise A0_ — (f>a — 4>b of the phase difference is then 
given by the number-phase Heisenberg inequality 



AiV_A(/)_ > 1. 



(3) 



In conventional interferometry, one of the inputs, say b, 
is a vacuum mode and the result of the interference is a 
binomial probability law, which yields the SNL, whatever 
the state of the bright input, a, is 

In order to achieve HLI, it is necessary to modify 
the vacuum input b. This can be done by squeezing 
a vacuum-field quadrature, as first proposed by Caves 
y and experimentally realized by Xiao et al. and 
Grangier et al. Ja. Numerous other schemes have been 
proposed since jj, in particular by Yurke, McCall, and 
Klauder and Holland and Burnett Q • 

This paper focuses on the latter proposal, whose ex- 
perimental implementation is in progress in our group 
Hi 0,0]. It is based, for reaching the HL, on the use of 
indistinguishable twin Fock states 



\ka,uj,e;n) ig) \kb,Lu,c,n), 
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where k is the unit wave vector, uj the frequency, e the 
unit polarization vector, and n the photon (boson) num- 
ber. It is remarkable that, in the same manner as for the 
SNL , common- mode quantum statistics do not matter 
and the HL can be reached as efficiently by using super- 
positions or mixtures of twin Fock states p^ . or twin 
beams, whose density matrix is of the form 



FIG. 1: Beam splitter. Reflection/transmission coefficients: 
r = —it = 2~^^^. a,b,c,d: annihilation operators. The wave 
fronts overlap perfectly, i.e. the beams are aligned so that 
their wave vectors verify fcc,d ~ ki,^a. 
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with all quantum numbers identical between modes a 
and &, except the wave vectors — which will be exactly 
matched at the beam splitter (Fig.QJ. 

Equations H4I5() describe perfectly correlated interfer- 
ometer inputs. The physics of the HL in this case lies 
solely in the generalized Hong-Ou-Mandel interference 
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[T^ of the twin Fock states Q at a beam splitter. Su- 
perpositions or mixtures © do not alter the result [l^ . 

This has been demonstrated experimentally using 
a pair of independently number-squeezed picosecond- 
pulsed beams from optical fibers jl4| and, more recently, 
by our group with ultrastable milliwatt-level twin beams 
emitted by a type II OPO electronically phase-locked at 
frequency degeneracy ^| . 

Good approximations of twin boson modes have also 
been implemented using internal atomic states of trapped 
ions |il5i | and other proposals include internal atomic 
states of Bose-Einstein condensates ITgI [iTl and orthog- 
onal phonons in a single trapped ion [l^- In the follow- 
ing, we will consider the case of photons, but the reader 
should bear in mind that the generalization to bosons 
is mathematically immediate. Some physical analogs of 
the optical beam splitter are well known: for example, 
a 7r/2 resonant laser pulse coupling two atomic energy 
states is isomorphic to a balanced beam splitter overlap- 
ping two optical paths. Both are SU(2) systems and can 
be described by the Schwinger representation [T^. 

A particularity of the Holland-Burnett scheme lies in 
its nonconventional detection process. When a twin 
mode input is used, the expectation values for the out- 
put photon numbers of an interferometer are independent 
of the phase difference </)_ H, 0, E| . Direct detection 
of the average output intensities, i.e. of the second-order 
moment of the input field, cannot therefore be used. The 
fourth-order moment is phase dependent, as has been 
noted first by Bouyer and Kasevich jl7j |. however, a sim- 
ple derivation of the variance of that measurement shows 
that its signal-to- noise ratio is very low: V2 at best , 
which lowers detection accuracy. 

Coincidence detection can dramatically improve the 
signal-to-noise ratio of a fourth-order moment measure- 
ment and a proof-of-principlc experimental demonstra- 
tion has been carried out for iV = 2 It is, how- 
ever, not clear how to extend coincidence measurements 
to photon numbers N > 2. What's more, an analysis 
of coincidence measurements in the presence of detector 
losses has shown extremely unfavorable results: in par- 
ticular, it seems that the losses need to be smaller than 
for this approach to yield even sub-SNL results [2ll| . 
These unfavorable results have cast a doubt on the ex- 
perimental feasibility of twin- mode HLI . 

However, we would like to establish here that such 
is not the case if one considers the original measure- 
ment strategy proposed by Holland and Burnett, which 
does not rely on coincidence detection but on a Bayesian 
reconstruction procedure from dynamical photocurrent 
measurements Q. The present paper thus focuses on 
Bayesian twin-mode interferometry in the presence of de- 
tection losses, and on its scaling with the photon number. 

The influence of losses on Bayesian HLI was investi- 
gated by numerical simulation in a previous paper [T^. 
but the results were limited to low photon numbers (100) 
and did not address the question of the scaling of the 
sensitivity with N. In this paper, we carry out numerical 



simulations of Bayesian measurements for photon num- 
bers ranging from iV = 10 to 10 000, in the presence of 
arbitrary detector losses. The results are very different 
from the ones obtained for coincidence detection and the 
prospects for experimental realizations are not as unfa- 
vorable as one might have feared from the conclusions of 

El- . 

In Section II of this paper, we briefly review quantum 
interferometry and the Bayesian measurement procedure. 
We then present our numerical simulations in Section III 
and conclude. 



II. TWIN-MODE HLI 

All optical interferometers can be modeled using a 
Mach-Zehnder model (Fig. O, i.e. a first beam splitter 




FIG. 2: Mach-Zehnder interferometer. 

"splits" the input light beam into two paths a and b 
which experience a phase difference to be measured and 
the fields are then "recombined" by a second beam split- 
ter. The goal is to measure this phase difference with 
high precision. As mentioned earlier, the output statis- 
tics of the first beam splitter determine the sensitivity 
and noise properties of the interferometer. 

For conventional interferometry, the input state 

\ka,u;,e;N) (g)\k,0J,e;0), (6) 

yields a binomial law for the beam splitter output statis- 
tics AA^_ oc a/^V If the input is a minimum uncer- 
tainty state — equality in — the standard deviation 
of the phase difference is then the SNL (|2Jl. 

A squeezed-vacuum input state can be used too: 

\ka,UJ,e;a) (g)\kb,uj,e;0,r), (7) 

where r is the squeezing parameter. If the squeezing is in 
quadrature with the coherent state phase-difference 
squeezing is obtained at the output of the beam splitter 

We now turn to twin-mode Heisenberg-limited inter- 
ferometry, where the input state is that of Eq. Q . The 
detail of the theory is given in and briefly recalled 

here. As mentioned in the Introduction, the expectation 
value of the output intensity of a Mach-Zehnder interfer- 
ometer is independent of the phase shift </> between the 
two arms. The Bayesian detection method is therefore 
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used. It consists of burst measurements of the output in- 
tensities, rather than average measurements yielding the 
expectation value. The probabili ty d istribution of the 
output intensity is easy to derive |l2| in the Schwinger 
representation |19j and is given by 

P{j,m,cl>)^^J—^[Pp{cos4>)f, (8) 
[j + m)] ■' 

where j — ijia + nfj)/2 and m = {ua — ni,)/2, Ua.b are the 
measured photon numbers at each output, and where 
PJ"(cos0) is an associated Legendre polynomial. Note 
that, if there are no losses, j = N always for input state 
Q because of unitarity. 

Assuming each measurement result independent of one 
another, the probability of a set of k burst measurements 
[(ji, TOi); . . . ; {jk,mk)] is simply the product 



(9) 



When plotted against (j) (Fig. |3Jl, -P'*^-* is fairly broad at 
k — 1 but narrows down quickly to a sharp peak with 
~ at k = 10. In the presence of optical losses, 
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FIG. 3: Simulated Bayesian twin-mode phase probability dis- 
tribution for a phase of 0.0001 radians (denoted by the verti- 
cal line) and several burst measurements ("modes") of 10000 
photons each, at 100% efficiency. 

analytical treatment is very complicated and numerical 
simulation is necessary, in particular in order to elucidate 
the scaling of the phase error with the photon number. 
This is the subject of the next section. 



the program samples a number of one-shot instantaneous 
intensity difference measurements at the photodetectors. 
This corresponds to a number difference m between the 
output beams. While the average number difference is 
expected to be zero over long measurement times, instan- 
taneous fluctuations are not. Thus in any given measure- 
ment the number difference may be drastically different 
from zero. Using this difference the program reconstructs 
a probability distribution for the phase, for each burst 
measurement of the number difference m. 

Each measurement gives a different m, most of the 
time, and can also give a different j if detectors are not 
100% efficient or the input state is not perfectly corre- 
lated. Previous work showed that both cases are equiv- 
alent, even though the nonideal-detector case is much 
easier to handle numerically, using Legendre in lieu of 
general Jacobi polynomials [13 ■ We therefore claim that 
a general analysis can be carried out by solely account- 
ing for detection losses and assuming ideal initial corre- 
lations. 

The "measured" peak of Fig. |31 is centered on the 
"true" phase value, fixed by the program, and the peak 
width gives the measurement uncertainty. The code cal- 
culates this uncertainty and shows a diminishing value 
as the number of bursts in the measurement increases. 
This occurs even if the peak predicted by the proba- 
bility distribution is not in fact centered at the proper 
phase value (postulated by the program before measur- 
ing the intensity difference). In such a situation, the 
distance of the measured peak from the real phase value 
becomes the dominant factor in the uncertainty, as the 
width of the distribution is generally much smaller than 
this value. On average this distance is much smaller than 
the SNL for all types of quantum efhciencies studied, and 
it reaches the HL for 100% efficiency. 
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III. NUMERICAL SIMULATIONS OF 
NONIDEAL BAYESIAN TWIN-MODE HLI 

A. Method 

The monte carlo simulation for this study uses twin 
Fock states as the input. The detectors along with their 
efficiencies are simulated on the output ports. For each 
phase difference measurement within a specified range. 



FIG. 4: Simulated Bayesian twin-mode phase probability dis- 
tribution for a phase of 0.01 radians (denoted by the vertical 
line) and several burst measurements ("modes") of 100 pho- 
tons each, at 99% efficiency. 

Several total input photon numbers were used in the 
program, 10, 100, 1000, and 10000 at efficiencies of 90%, 
99%, 99.9%, 99.99%, and 100%. 

For each phase difference specified there is an associ- 
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ated Legendre polynomial evaluated that is used to simu- 
late the detected total photon number and number differ- 
ence between the detectors. The losses are computed by 
operating on the total number of photons in each beam as 
they exit the second beam splitter using a binomial dis- 
tribution weighted by the detector efficiency. The total 
number of photons on the output ports is simulated using 
a random possibility to transmit or reflect each photon 
at the beam splitters, which depends on the phase. 

For examples of the probability distributions, or syn- 
onymously, phase distributions, see Figs. El and 01 (note 
the vertical axes are in arbitrary units that indicate rel- 
ative probability compared to points on the same phase 
distribution). 
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FIG. 5: Average of measured minus real phase for 100% effi- 
ciency (log-log scale). 



B. Results 

All of the figures in this section refer to uncertainties 
on measurements of phase angles (j) = 1/N. This means 
the data points in each graph refer to the uncertainty 
in a measurement for a different phase angle from one 
data point to the next. The most important results are 
found for phase differences ranging from zero to a few 
times 1/N. There is no "optimal" phase angle in this 
range, unlike what Kim et al. found for coincidence detec- 
tion |23|. For large phase angles {(j) 7r/2), one cannot 
acheive better than the SNL. Further, the uncertainty is 
found to be quasi-ideal for nonideal quantum detection 
efficiencies = rjN such that 



(10) 



We have also found the remarkable result that the phase 
measurement uncertainty scales as 1 /N for such nonideal 
quantum efficiencies, as in the case of 100% efficient de- 
tectors. This defines a margin of experimental conditions 
within which the phase uncertainty may be further in- 
creased below the SNL by merely increasing the boson 
number N, but not the efficiency itself (or the effective 
photon correlation). Even though rjN may be thought of 
impossible to reach at optical powers, such may not be 
the case with matter waves, where very high efficiency 
detection is possible. 

While detector efficiencies smaller than 77 at do not yield 
such a favorable scaling of the phase measurement er- 
ror (1/VlV instead of 1/iV, like squeezing measurements 
|23|), it is nevertheless found that the Bayesian measure- 
ment method can always give a better precision than the 
SNL, which, again, was not the case for coincidence de- 
tection |2lj |. In this regime, the gain in phase error is 
solely conditioned by the improvement of the effective 
correlation, as is the case for vacuum squeezing experi- 
ments p^ . 

We now turn to the simulation plots. For 100% effi- 
ciency (Fig. EJ the phase error follows 1/A^ as expected 
and one can also see that both the simulated SNL and 
HL are slightly below or 1/A^. We attribute this 



to the averaging process of multiple-burst measurements. 
In situations involving linear averaging, one gets a reduc- 
tion factor of the noise of 1/ y/p, where p is the number 
of data samples. In our case, the measurement probabil- 
ities are multiplied together, not averaged. Nevertheless, 
the results are found to lay within less than an order of 
magnitude of a line consistent with the HL reduced by 
1 / where p is the number of probability distributions, 
or bursts. 

For nonideal efficiencies 77 < 100% (Figs. to |SJ), we 
find that the phase error follows 1/A^ for low photon num- 
bers before crossing the HL line to follow the SNL at 
l/\/lV. Our simulations indicate that the inflexion point 
is 1/(1 - r/). For cases of 90%, 99%, and 99.9% effi- 




10 



100 1000 
photon number 



10000 



FIG. 6: Average of measured minus real phase for 99.9% ef- 
ficiency (log- log scale). 

ciency, as one would expect, scaling occurs as 1/A^ for 
small j, and then rapidly depreciates as photon number 
increases. This can be explained by the probability dis- 
tribution being dependant on the detector efficiency, i.e. 
the losses in the total photon number. At 99% for 100 
photons one would expect to detect on average 99 pho- 
tons per measurement, while measuring 10000 photons at 
an efficiency of 99% naturally produces much larger gross 
losses, even though the percentage losses are equivalent 
for the two cases. The case of 90% efficiency is shown in 
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FIG. 7: Average of measured minus real phase for 99% effi- 
ciency (log- log scale). 



FIG. 9: Average of measured minus real phase for 1 — 1/A^ 
efficiency. 
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C. Interpretation 



As the phase angle measured here suggests, in typical 
experiments the measured phase will need to be within 
a small range. As phase angle dilference increases one 
obtains larger error bars on the measured phase, as fig- 
ure^] indicates. This is true for all efRciencies studied 
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FIG. 8: Average of measured minus real phase for 90% effi- 
ciency. 



FiglHlto illustrate this point. Important to note is that 
the average error stays below the SNL, even for 10000 
photons. Between 500 and 10000 one can observe the 
scaling take on a definite nature of 1/ \/N. The distance 
below the classical line is within an order of magnitude 
and is again due to the use of Bayes theorem, coupled to 
the degree of effective correlation determined by the loss 
level. It is also important to note that this line is well 
below the "true" SNL, which is given by 
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FIG. 10: Average of measured phase vs real phase for 100 
photons at 99% efficiency. 



All other efficiencies naturally follow this trend. 

This leads us to finally consider the particular quantum 
efficiency rjN oi Eq. ^| The results of the simulations 
are plotted in Fig. |21 which shows, indeed, a behavior 
comparable to that of the ideal efficiency of Fig. O and 
has been confirmed for phases within the range of = 
to several times (j) = 1/N. 



here. As is clearly evident from the graph, the error bars 
at higher phase differences are much greater in magni- 
tude than the distance of the measured phase from the 
line. The opposite case prevails at low phase differences. 
Kim et al have suggested that an optimal phase angle ex- 
ists for each efficiency for which uncertainty is minimized 
PH. However, if one opts for the Bayesian postprocess- 
ing of data analyzed here, rather than for coincidence 
measurements, then evidently there is no optimal phase 
difference, and any minute phase shift around a multiple 
of TT radians can be measured with high confidence. 
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IV. CONCLUSION 

We have presented a theoretical study of twin Fock 
states incident on a Mach-Zehnder interferometer, with 
simulated measurements based on the Bayesian probabil- 
ity reconstruction first proposed by Holland and Burnett 
m 01 . The lingering question of the scaling of the phase 
measurement error with the photon number has been elu- 
cidated and it has been found that: 

• One may beat the SNL for nonideal detector effi- 
ciencies without the signal-to-noise impediment of direct 
detection of the fourth-order correlation 0, 0| . This 
paves the way to the realization of a sub-SNL interfer- 
ometry experiment with phase- locked twin beams 
an interesting prospect since intensity-difference squeez- 



ing is one of the most performant experimental systems 
available. 

• Nonideal efficiencies can still allow one to reach the 
HL provided that the losses are of the order of 1 /N . Such 
requirements may be too stringent for photon interfer- 
ometry, but Bose-Einstein condensates or other atom or 
trapped-ion experiments may be good candidates. 

• The problems found with coincidence detection [2ll | 
are not relevant to Bayesian twin-mode HLI. 
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